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ABSTRACT 

We have carried out a multi-site campaign to measure oscillations in the F5 star Procyon A. We 
obtained high-precision velocity observations over more than three weeks with eleven telescopes, with 
almost continuous coverage for the central ten days. This represents the most extensive campaign so 
far organized on any solar-type oscillator. We describe in detail the methods we used for processing 
and combining the data. These involved calculating weights for the velocity time series from the 
measurement uncertainties and adjusting them in order to minimize the noise level of the combined 
data. The time series of velocities for Procyon shows the clear signature of oscillations, with a plateau 
of excess power that is centred at 0.9 mHz and is broader than has been seen for other stars. The mean 
amplitude of the radial modes is 38.1 ± l.Scms" 1 (2.0 times solar), which is consistent with previous 
detections from the ground and by the WIRE spacecraft, and also with the upper limit set by the 
MOST spacecraft. The variation of the amplitude during the observing campaign allows us to estimate 
the mode lifetime to be l.S^g'g d- We also find a slow variation in the radial velocity of Procyon, with 
good agreement between different telescopes. These variations are remarkably similar to those seen 
in the Sun, and we interpret them as being due to rotational modulation from active regions on the 
stellar surface. The variations appear to have a period of about 10 days, which presumably equals the 
stellar rotation period or, perhaps, half of it. The amount of power in these slow variations indicates 
that the fractional area of Procyon covered by active regions is slightly higher than for the Sun. 
Subject headings: stars: individual (Procyon A) — stars: oscillations 
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1. INTRODUCTION 

Measuring solar-like oscillations in main-sequence and 
subgiant stars requires high-precision observations - ei- 
ther with spectroscopy or photometry - combined with 
coverage that is as continuous as possible. Most of 
the results have come from high-precision Doppler mea- 
surements using ground-based spectrographs, while mea- 
surements from spacecraf t h ave also been repo rted (see 
iBedding fe Kieldsenll2007l and lAerts et al.ll2008l for recent 
summaries). 

Procyon has long been a favourite target for oscil- 
lation searches. At least eight separate velocity stud- 
ies have reported an e xcess in the power s pectrum, be- 
ginning with that by iBrown et al.l (|199lD . which was 
the first report of a solar-like power exce ss in another 
star. For the most recen t exam p les, see | Martic et al.l 
(12004 . lEggenberger et aT] (|2004l ). iBouchv et all (|2004f ) 
and IlLeccia et al.l (|2007f h These studies agreed on the 
location of the excess power (around 0.5-1.5 mHz) but 
they disagreed on the individual oscillation frequencies. 
However, a consensus has emerged that the large sep- 
aration (the frequency separation between consecutive 
overtone modes of a given angular degree) is about 
55 /iHz. Ev i dence for this value was first given by 
iMosser et al.l (Il998f) and th e first clear detection was 
made bv lMartic et al.l (| 19991 ). 

Controversy was generated when photometric obser- 
vations obtained with the MO ST satellite failed to 
reveal evidence for o s cillations dMatthews et alJ l2004t 
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20071: iBaudin et alJ I2008D . However 



20051 ) argued that the MOST non- 



detection was consistent with the grou nd-based data. 
Meanwhile. iRegulo fc Roca Cortes! (( 20051) suggested that 
the signature of oscillations is inde ed present in th e 
MOST data at a low level (see also lMarchenkdl2008f ). 
Using space- based p hotometry with the WIRE satellite, 
iBruntt et al.l (|2005f) extracted parameters for the stellar 
granulation and found evidence for an excess due to os- 
cillations. 

All published velocity observations of Procyon have 
been made from a sing le site, with the ex ception of 
two-site observations by iMartic et alJ f2004) . Here we 
describe a multi-site campaign on Procyon carried out 
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in 2007 January, which was the most extensive velocity 
campaign so far organized on any solar-type oscillator. 
The only other comparable effort to measure oscillations 
in this type of star was the mu lti-site photometric c am- 
paign of the open cluster M67 (|Gilliland et al.lll993l ). 

2. VELOCITY OBSERVATIONS 

We observed Procyon from 2006 December 28 until 
2007 January 23, using a total of eleven telescopes at 
eight observatories. These are listed in Table [TJ ordered 
westward by longitude. Note that the FIES spectrograph 
on the Nordic Optical Telescope was still being commis- 
sioned during the observations and the velocity precision 
is therefore somewhat lower than for the other telescopes. 

The team members from each telescope were respon- 
sible for producing a velocity time scries from the ob- 
servations, together with estimates of uncertainties. In 
six of the spectrographs, the stellar light was passed 
through an iodine absorption cell to provide a stable 
wavelength reference. In four others, wavelength cali- 
bration was achieved by recording the spectrum from a 
thorium-argon emission lamp alongside the stellar spec- 
trum, while with EMILIE, exposures of the stellar spec- 
trum were alternated with those of a white-light source 
passing through an iodine cell. Details of the methods 
used with each spectrograph are given in the references 
listed in Table [1] and details of the observations are given 
in Tabled 

The velocity time series of Procyon is shown in Fig.[T]a, 
using a different color for each telescope. Differences be- 
tween telescopes in the absolute zero point of velocity 
are not significant, and so all the curves have been shifted 
into alignment by subtracting a constant offset. This was 
done by setting the velocities from each telescope to have 
zero mean, excepting Lick and SARG, for which better 
alignment was achieved by setting the means to —6 m s -1 
and -5 ms -1 , respectively. Note that EMILIE is not 
shown in Fig. [TJa because those data were referenced to 
a different value on each night and so the night-to-night 
variations are not measurable (this does not affect their 
usefulness for oscillations studies, however). FIES is also 
not shown in Fig.[T]a because it has much greater scatter 
than the rest. 

In Fig. [TJa we see variations in the radial velocity of 
Procyon on timescales of days. The good agreement be- 
tween the different telescopes indicates that these slow 
variations have a stellar origin, although the imperfect 
match in overlapping sections shows that there is also a 
contribution from instrumental drifts. Figure [If) shows 
a close-up of the central part of the campaign, during 
which the coverage was above 90%. The solid curve 
shows the velocities after smoothing, to better reveal the 
slow variations, which are discussed in ij4. II below. 

While interesting in their own right, the slow variations 
in the velocity series significantly affect our ability to de- 
tect oscillations, due to spectral leakage of power from 
the low-frequency part of the spectrum to the oscillation 
region at higher frequencies. Before merging the data 
from the individual telescopes, we therefore removed the 
low-frequency variations from the velocity series. This 
was done for each telescope by removing all the power 
below 280/iHz, a value that was chosen so as to effec- 
tively remove the slow variations without affecting the 
oscillation signal. This filtering was done by subtracting 
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a smoothed version of the time series that contained all 
the power below that cut-off frequency. 

In Fig[T]c we show a close-up of a segment during which 
three spectrographs were observing simultaneously. The 
stellar oscillations are clearly visible, with typical peri- 
ods of about 15 minutes, and there is good agreement 
between the different telescopes. Note that these data 
have been filtered to remove the slow variations. 

3. OPTIMIZING THE WEIGHTS 

The procedures for extracting velocities for each tele- 
scope also produced estimates of the uncertainties, Ci. 
In our analysis, we used these uncertainties to calcu- 
late noise-optimized weights in the usual way, namely 
Wi = l/of . If weights are not used when calculating the 
power spectrum, a few bad data points can dominate and 
increase the noise floor significantly. 

We now describe the process we used to adjust these 
weights, which aims to minimize the noise level in the 
final power spectrum. The procedure involves identify- 
ing and revising those uncertainties that were too op- 
timistic, and at the same time rescaling the uncertain- 
ties to be in agreement with the actual noise levels in 
the data. These met hods have already been described 
m previous paper s (IButler et all l2004t iBedding et all 
l2007t iLeccia et al.ll2007ft . but the present analysis differs 
slightly from those descriptions and we will therefore de- 
scribe them in some detail. One difference is that the 
analysis had to be tailored to the individual time series 
because of the large range of Nyquist frequencies (see 
Table p. 

To illustrate the process, we show in Figs. [2] and [3] 
segments of data at different stages in the process for 
two telescopes (HARPS and EMILIE). The top panels 
(Figs. EJa and [St) show the velocities for a single night, 
with the slow variations removed. The remaining panels 
show the uncertainties at different stages in the analysis. 

It is important to stress that we are not adjusting the 
velocities, only the uncertainties. Of course, those ad- 
justments still affect the power spectrum of the velocities 
(which is, after all, why we are making the adjustments) 
and so it is important to ensure that they do not distort 
the oscillation signal and that the final weights reflect as 
accurately as possible the actual noise properties of the 
series. 

3.1. Scaling the uncertainties 

We have scaled the uncertainties so that they agree 
with the noise level in the corresponding amplitude spec- 
trum, (Tamp, as measured at high frequencies. This was 
done for each night and each telescope by multiplying 
the uncertainti es, gj, by a constant so that they satisfied 
equation (3) of IButler et all (|2004l ): 

N 

Camp ^ = W 
i=l 

This scaling was repeated after each step in the process 
described below. Figures and[3t show the uncertain- 
ties after scaling, and before any further adjustments. 

3.2. Filtering the uncertainties 

It is clear that the uncertainties in some parts of the 
time series show variations that correlate with the oscil- 
lations of Procyon. The clearest example is HARPS, as 



shown by comparing the top two panels of Fig. [2l but 
the effect is also visible for other telescopes (e.g., Fig[3|). 
To remove this structure in the uncertainties, we have 
bandpass-filtered each of the uncertainty time series to 
remove all power in the frequency range 280-2200/iHz. 
This removed fluctuations in the weights on the timescale 
of the stellar oscillations, while retaining information on 
longer timescales (such as poorer conditions at the begin- 
nings and ends of nights) and on shorter timescales (such 
as individual bad data points). This process resulted in 
slightly lower noise levels in the final power spectra for 
some of the individual telescopes, reflecting the fact that 
the uncertainties, like any measurement, contain noise 
that is reduced by bandpass filtering. Figures HJc and[3]s 
show the uncertainties after filtering. 

We also noticed a few data points (~10) with unreal- 
istically low uncertainties. These points would be given 
too high a weight in the analysis and the uncertainties 
were therefore reset to the mean uncertainty for that tele- 
scope night. 

3.3. Down-weighting bad data points 

Down- weighting of bad dat a points was done follow- 
ing the method described by IButler et al.l ((2004) , with 
one difference that is discussed below. The first step was 
to make a high-pass-filtered version of the velocity time 
series in which both the slow variations and the stellar 
oscillations were removed. This gave us a series of resid- 
ual velocities, fj, in which we could identify data points 
that needed to be down-weighted, without being affected 
by spectral leakage from the oscillations. The frequency 
limit of this high-pass filter varied from telescope to tele- 
scope, depending on the Nyquist frequency of the data. 

We compared the velocity residuals, fj, with the cor- 
responding uncertainty estimates, er^. Bad data points 
are those for which the ratio \rijai\ is large, i.e., where 
the residual velocity deviates from zero b y more than 
expect ed from the uncertainty estimate. IButler et al.l 
(2004), who analyzed data similar to ours, found that 
the fraction of good data points was essentially unity up 
to \ri/<Ji\ — 2 and then dropp ed off quickly for lar ger val- 
ues of \n/o-i\ - see Figure 3 of IButler et "all (|2004h . They 
therefore introduced the factor /, which is the fraction of 
good data points as a function of | r f / cr^ | and which they 
obtained as the ratio between the distribution of data 
points in a cumulative histogram of \ri/o~i\ and a best-fit 
Gaussian distribution. 

We used a slightly different approach. With the knowl- 
edge that points with large values of |r»/<7i| are bad, we 
introduced an analytical function 

f( x i) = r x 10 ' Xi = ktMli ( 2 ) 



which has shape v ery similar to the fr action / as a func- 
tion of |rj/<Tj| in IButler et all (|2004l ). The adjustable 
parameter xq controls the amount of down-weighting; it 
sets the value of \ri/ai\ for which the weights are multi- 
plied by 0.5, and so it determines how bad a data point 
should be before it is down-weighted. The optimum 
choice for xo was found through iteration, as described 
below. Once this was done, we used f(xi) to adjust 
the we ights by dividing er,; by y/ f(xi), as in IButler et al.l 
(2003). 
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The noise level used for optimizing xq was measured 
in a frequency band near 2 mHz in a weighted ampli- 
tude spectrum, between the oscillations and the high- 
frequency part of the spectrum used for determining the 
|r»/<7{| values. The exact position of the band was cho- 
sen for each spectrograph separately, because of the dif- 
ferences in Nyquist frequencies. For each trial value of 
Xo, the noise level was determined from a time series in 
which all power had been removed at both the low- and 
high-frequency side of the frequency band, i.e., from a 
bandpass-filtered time series containing information in 
the specific frequency band only. This was done because 
applying the weights changes the spectral window func- 
tion and thus the amount of spectral leakage into the 
frequency band where we determine the noise: if we did 
not filter out the oscillations and the high-frequency part 
of the spectrum, we would optimize for a combination of 
low noise and minimum amount of spectral leakage (from 
both the low- and high frequency side of the passband). 
In other words, the spectral window function would influ- 
ence our choice of xq, which is not optimal for obtaining 
the lowest possible noise level. 

The procedure described above was repeated for a 
range of xq values and we chose the one that resulted 
in the lowest noise in the power spectrum. Depending 
on the telescope, and hence the noise properties of the 
time series, the optimal values of xq ranged from 1.7 to 
4.3. Figures and[3]ff show the final uncertainties for 
HARPS and EMILIE. 

This completes our description of the process used to 
adjust the uncertainties. The results from calculating 
weighted power spectra using these uncertainties are pre- 
sented in §4.21 First, however, we discuss the slow vari- 
ations in the velocity of Procyon that are present in our 
data. 

4. RESULTS 

4.1. Slow variations in stellar velocity 

The slow variations in the radial velocity of Procyon 
seen in Fig. [1] are remarkably similar to those seen in 
the Sun. This can be seen from Fig. [H which shows a 
typical time series of solar velocity measurements made 
with the GOLF instrument on the S OHO spacecraft 
(lUlrich et all l2000t IGarcia et all l2005h . Lon ger series 
of GOLF data show variations with a period of about 
13 d arising from active region s crossing the solar disk 
(Fig 11 in IGarcia et al.l 120051: see also iTinglev et all 
2008). This 13- day periodicity in so lar velocities was 
first observed by IClaverie et al.1 (|1982t ). who attributed 
it to rapid rotation of the core, but it was subse- 
quently shown to be due to surface rotation of active 
regions dDurrant k. Schroterl 119831 : lAndersen &i Maltbvl 
fl98l lEdmunds fc Gougr]ll983f) ~imilarlv. we attribute 
the slow variations in the radial velocity of Procyon to 
the appearance and disappearance of active regions and 
their rotation acro ss the stellar disk. This explanation 
was also invoked bv lMosser et al.l (|2005f ) to explain much 
larger velocity variations measured with HARPS for the 
COROT target HR 2530 (= HD 49933; spectral type 
F5V). 

The slow variations in Procyon appear to have a pe- 
riod of Psiow = 10.3±0.5d, as measured from the highest 
peak in the power spectrum. This agrees with an appar- 
ent periodicity of about 10 d in two-site observations of 



Procyon obtained over 20 nights bv lKambe et~aT1 (|2008l . 
see their Fig. 12). Identifying this as t he stellar rota- 
tion p eriod and using a radius of 2.05 P© (jKervella et al.1 
1200 if ) implies a surface rotational speed at the equator 
of v = 10 ± 0.5 kms -1 . The measured value of P s iow 
perio d might also correspond to half the rotation pe- 
riod (|Clarkd [20031. in which case the rotational speed 
would be half the value given above. Combining with 
the spectroscopic value of v sin i = 3.16 ± 0.50kms _1 
(|Allende Prieto et al.l 12002) gives an inclination angle of 
i = 18 ± 3° (if P rot = P s iow) and i = 39 ± 7° (if 
P ro t = 2P s i ow ). The inclinat ion of the binary orbit is 
31.1 ± 0.6° (jGirard et alJl2000f ) and so, if we require that 
the rotation axis of Procyon is aligned with the orbital 
rotation axis, it may be that P rot = 2P s i ow . 

4.2. Power spectra and stellar activity 

The weighted power spectrum for each telescope and 
for the combined time series, based on the uncertainties 
discussed in §21 are shown in Fig. [5] The noise levels, as 
measured at high frequencies in the amplitude spectrum 
(camp)i are given in column 8 of Table [2] The final col- 
umn of that table gives the mean noise level per minute 
of observing time, with a spread that reflects a combi- 
nation of factors, including telescope aperture, observing 
duty cycle, spectrograph design, as well as atmospheric 
conditions such as seeing. 

The power spectrum of the combined time series is 
shown again in Fig. [6l both with and without the use of 
weights. As is well-established, using weights reduces the 
noise level significantly, at the cost of increasing the side- 
lobes in the spectral window (because the best data seg- 
ments are given more weight - see insets) . When weights 
are used, the noise level above 3 mHz is 1.9 cm s -1 in 
amplitude, but this does include some degree of spectral 
leakage from the oscillations. If we high-pass filter the 
spectrum up to 3 mHz, the noise level drops to 1.5 cm s _1 
in amplitude. Note that without the use of weights, the 
noise level is higher by more than a factor of two. 

Looking again at Fig. [6l we see that the use of weights 
appears to have increased the amplitude of the oscilla- 
tions. In fact, this indicates the finite lifetime of the 
oscillation modes: in Fig. [6^ the HARPS data are given 
the highest weight, and so the effective duration of the 
observations is decreased (and the sidelobes in the spec- 
tral window become much stronger). Our estimate of the 
mode lifetime is given in £14.31 

It is also useful to convert to power density, which is in- 
dependent of the observing window and therefore allows 
us to compare noise levels. This is done by multiplying 
the power spectrum by the effective length of the observ- 
ing run, which we calculate as the reciprocal of the area 
integrated under the spectral window (in power). The 
values for the different telescopes are given in column 3 
of Tabled In Fig. [7] we show the power density spectrum 
on a logarithmic scale for the HARPS data, which has the 
lowest noise per minute of observing time. We see three 
components: (i) the oscillations (about 300-1 100 /iHz); 
(ii) white noise at high frequencies; and (iii) a sloping 
background of power at low frequencies (stellar granula- 
tion and activity, and presumably also some instrumen- 
tal drift). Figure [8] compares the power density spectra 
for the different telescopes. They show a similar oscilla- 
tion signal and similar background from stellar noise at 
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lower frequencies (below about 250 /xHz), with different 
levels of white noise at higher frequencies (above about 
2000 /iHz), reflecting the different levels of photon noise. 

In Fig. the lower two dashed lines indicate the back- 
ground level in the Sun, as measured from the GOLF 
data during solar minimum and maximum, respectively. 
The upper dashed line is the solar maximum line shifted 
to match the power density of Procyon, which required 
multiplying by a factor of 40. We can use this scal- 
ing factor to estimate the fraction of Procyon's surface 
that is covered by active regions, relative to the Sun, as 
follows. The low-frequency part of the velocity power- 
density spectrum from the Sun falls off as frequency 
squared (|Harvevlll985HPalle et a l. 1999), and we see the 
same behaviour for Procyon. Hence, in both cases we 
have 

PD{v) oc v- 2 . (3) 

Let T be the typical time for an active region to be vis- 
ible on the surface (which may depend on both rotation 
and the typical lifetime of active regions). The ampli- 
tude of the velocity signal at frequency vq = l/T will 
be proportional to the fractional area covered by active 
regions, da/a, and to the projected rotational velocity, 
vsini. The power density at vq is therefore: 

PD(y Q ) = (da / a) 2 (v sin i) 2 . (4) 

Combining these gives 



, da 
PD{ V ) = — 
1 a 



v sin i 
T 



(5) 



We will assume that T is proportional to the rotation pe- 
riod, which we take to be 10.3 d (or twice that value) for 
Procyon and 25.4 d for the Sun. The measured values for 
v sin i are 3.2 km s -1 for Procyon (jAllende Prieto et alJ 
2002) and 2.0kms _1 for the Sun. Combining these val- 
ues with our measurement of the power densities indi- 
cates that the area covered by active regions on Procyon 
is about 1.6 times the solar maximum value (or twice 
that value). No detection of a magnetic field in Procyon 
has been reported, and published upper limits imply that 
the average fiel d cannot be more than a few times solar 
(see Table 3 in iKim et alJ 12007). which appears to be 
consistent with our results. 

4.3. Oscillation amplitude and mode lifetime 

To measure the amplitude of oscillatio ns in Procyon, 
we ha ve used the method described by iKieldsen et al.l 
(2008). In brief, this involves the following steps: 

(i) heavily smoothing the power spectrum (by convolving 
with a Gaussian having a full width at half maximum of 
4Ai/, where Av is the large frequency separation), to pro- 
duce a single hump of excess power that is insensitive to 
the fact that the oscillation spectrum has discrete peaks; 

(ii) converting to power density (see ^4.2[) ; (iii) fitting 
and subtracting the background noise; and (iv) multi- 
plying by A^/4.09 and taking the square root, in order 
to convert to amplitude per radial oscillation mode. Note 
that 4.09 is the effective number of modes per order for 
full-disk velocities observations, no rmalized to the ampli- 
tudes of radial (I = 0) modes - see IKieldsen et all (2008) 
for details. 

We applied this method to each of the telescopes sep- 
arately, and the result is shown in Fig. [5J There are sig- 
nificant differences between the different curves, which 



we attribute to intrinsic variations in the star arising 
from the stochastic nature of the excitation and damping. 
To investigate this further, we also applied the method 
to the combined time series, after first subdividing it 
into ten 2-day subsets. Figure fTOl shows these amplitude 
curves and also their average. 

The amplitude curve of Procyon has a broad plateau, 
rather than the single peak that has been seen for 
other stars. Figure [IT] shows the smoothed amplitude 
curve for Procyon compared to the Sun and other stars. 
It is an updated version of Fig. 8 of IKieldsen et al.l 
(2008), where the followin g stars have been added: 
vl Ar a (iBouchv et alJ I2005T) HP 49933 (iMosser et al 
20051) u Her (iBonanno et al.ll2008h 7 Pav (jMosser et al 
2008bJ), r Cet (jTeixeira et al 112008ft 



The plateau for Procyon is centred at 900 /iHz and is 
about 500 /iHz wide, with a mean amplitude across that 
range of 38.1 ± 1.3 cms -1 . This is our estimate for the 
amplitude of the radial (I — 0) modes in Procyon. Com- 
paring with the co rresponding measure ment for the Sun 
(18.7±0.7cms -1 ; |Kjeldiin et al.ll2008[ ) implies that the 
velocity oscillations in Procyon are 2.04±0.10 times solar. 
In both Procyon and the Sun, the modes with I = 1 are 
highe r by a factor of 1.35 (see Table 1 of IKieldsen et al.l 
l2008h . 

The corresponding intensity amplitude, after account- 
ing for th e higher effective temperat ure of Procyon (see 
Eq. 5 in Ki eldsen fe Bedding] [1995D. is 1.60 times so- 
lar. This implies an amplitude at 500 nm of 6.8 ppm 
for I = and 8.5 ppm for I = 1 (see Table 1 of 
IKieldsen et al.ll2008l ). These amplitudes are completely 
consist ent with the detectio n of oscillations in Procyon by 
WIRE (jBruntt et a l. 2005) and with the upper limit s et 
by MOST dMatthews et al.ll200l iBedding et al.l l2005L 

The standard deviation of the ten segments in Fig.[Tt)lis 
a A / A = 10.4%±2.3%, which reflects the finite lifetime of 
the m odes. We can use equation (3) from IKieldsen et all 
(2008) to estimate the mode lifetime, but we must ac- 
count for the much greater width of the oscillation enve- 
lope in Procyon. Note that this equation was established 
empirically and we have confirmed it analyt ically using 
the work of iToutain fc Appourchauxl (|1994l) . We esti- 
mate the mode lifetime to be r = 1.5l '|d. This equals, 
withi n rather large uncer tainties, the solar value of 2.9 d 
(e.g., IChaplin et al.lll997l ). 

5. CONCLUSIONS 

We have presented multi-site velocity observations of 
Procyon that we obtained with eleven telescopes over 
more than three weeks. Combining data that spans a 
range of precisions and sampling rates presents a signifi- 
cant challenge. When calculating the power spectrum, it 
is important to use weights that are based on the mea- 
surement uncertainties, otherwise the result is dominated 
by the noisiest data. We have described in detail our 
methods for adjusting the weights in order to minimize 
the noise level in the final power spectrum. 

Our velocity measurements show the clear signature of 
oscillations. The power spectrum shows an excess in a 
plateau that is centred at 0.9 mHz and is broader than 
has been seen for other solar-type stars. The mean am- 
plitude of the radial modes is 38.lil.3cms- 1 (2.04±0.10 
times solar) , which is consistent with previous detections 
from the ground and by the WIRE spacecraft, and also 
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with the upper limit set by the MOST spacecraft. The 
variation of the amplitude during the observing cam- 
paign allowed us to estimate the mode lifetime to be 
1 c i ?1 - 9 (i 

J--O-0.8 Q - 

We also found a slow variation in the radial velocity 
of Procyon, with good agreement between different tele- 
scopes. These variations are remarkably similar to those 
seen in the Sun, and we interpret them as being due to 
rotational modulation from active regions on the stellar 
surface. The variations appear to have a period of about 
10 days, which presumably equals the stellar rotation pe- 
riod or, perhaps, half of it. The amount of power in these 
slow variations indicates that the fractional area of Pro- 
cyon covered by active regions is slightly higher than for 
the Sun. 

The excellent coverage of the observations and the high 
signal-to-noise should allow us to produce a good set of 



Aerts, C, Christensen-Dalsgaard, J., Cunha, M., & Kurtz, D. W., 

2008, Sol. Phys. in press ( arXiv : 0803 . 3527 ). 
Allende Prieto, C, Asplund, M., Lopez, R. J. G., & Lambert, D. L., 

2002, ApJ, 567, 544. 
Andersen, B. N., & Maltby, P., 1983, Nat, 302, 808. 
Baudin, P., Appourchaux, T., Boumier, P., Kuschnig, R., 

Leibacher, J. W., & Matthews, J. M., 2008, A&A, 478, 461. 
Bedding, T. R., & Kjeldsen, H., 2007, In Stancliffe, R. J., Houdek, 

G., Martin, R. G., & Tout, C. A., editors, Unsolved Problems 

in Stellar Physics: A Conference in Honour of Douglas Gough, 

volume 948, page 117. American Institute of Physics. 
Bedding, T. R., Kjeldsen, H., Arentoft, T., et al., 2007, ApJ, 663, 

1315. 

Bedding, T. R., Kjeldsen, H., Bouchy, F., et al., 2005, A&A, 432, 
L43. 

Bonanno, A., Benatti, S., Claudi, R., Desidera, S., Gratton, R., 
Leccia, S., & Paterno, L., 2008, ApJ, 676, 1248. 

Bouchy, F., Bazot, M., Santos, N. C, Vauclair, S., & Sosnowska, 
D., 2005, A&A, 440, 609. 

Bouchy, F., & Carrier, F., 2002, A&A, 390, 205. 

Bouchy, F., Maeder, A., Mayor, M., Megevand, D., Pepe, F., & 
Sosnowska, D., 2004, Nat, 432, 7015. 

Bouchy, F., Schmitt, J., Bertaux, J.-L., & Connes, P., 2002, In 
Aerts, C, Bedding, T. R., & Christensen-Dalsgaard, J., editors, 
IAU Colloqium 185: Radial and Nonradial Pulsations as Probes 
of Stellar Physics, volume 259, page 472. ASP Conf. Ser. 

Brown, T. M., Gilliland, R. L., Noyes, R. W., & Ramsey, L. W., 
1991, ApJ, 368, 599. 

Bruntt, H., Kjeldsen, H., Buzasi, D. L., & Bedding, T. R., 2005, 
ApJ, 633, 440. 

Butler, R. P., Bedding, T. R., Kjeldsen, H., et al., 2004, ApJ, 600, 
L75. 

Butler, R. P., Marcy, G. W., Williams, E., McCarthy, C, Dosanjh, 

P., & Vogt, S. S., 1996, PASP, 108, 500. 
Chaplin, W. J., Elsworth, Y., Isaak, G. R., McLeod, C. P., Miller, 

B. A., & New, R., 1997, MNRAS, 288, 623. 
Clarke, D., 2003, A&A, 407, 1029. 

Claudi, R. U., Bonanno, A., Leccia, S., Ventura, R., Desidera, S., 
Gratton, R., Cosentino, R., Paterno, L., & Endl, M., 2005, A&A, 
429, L17. 

Claverie, A., Isaak, G. R., McLeod, C. P., van der Raay, H. B., 
Palle, P. L., & Roca Cortes, T., 1982, Nat, 299, 704. 

Durrant, C. J., & Schroter, E. H., 1983, Nat, 301, 589. 

Edmunds, M. G, & Gough, D. O., 1983, Nat, 302, 810. 

Eggenberger, P., Carrier, F., Bouchy, F., & Blecha, A., 2004, A&A, 
422, 247. 

Endl, M., Cochran, W. D., Hatzes, A. P., & Wittenmyer, 
R. A., 2005, In Revista Mexicana de Astronomia y Astrofisica 
Conference Series, volume 23, page 64. 

Frandsen, S., & Lindberg, B., 2000, In Teixeira, T., & Bedding, 
T. R., editors, The Third MONS Workshop: Science Preparation 
and Target Selection, page 163. Aarhus: Aarhus Universitet. 

Garcia, R. A., Turck-Chieze, S., Boumier, P., Robillot, J. M., 
Bertello, L., et al., 2005, A&A, 442, 385. 



oscillation frequencies for Procyon. This analysis will be 
presented in subsequent papers. 



This work was supported financially by the Danish 
Natural Science Research Council, the Australian Re- 
search Council, the Swiss National Science Foundation, 
NSF grant AST-9988087 (RPB) and by SUN Microsys- 
tems. We thank Hugh Jones, Chris Tinney and the 
other members of the Anglo- Australian Planet Search 
for agreeing to a time swap that allowed our AAT ob- 
servations to be scheduled. MM is grateful to Prof. 
N. Kameswara Rao, G. Pandey and S. Sriram for their 
participation in campaign with VBT Echelle Spectrom- 
eter, which was used for the first time for Doppler spec- 
troscopy observations. 



Gilliland, R. L., Brown, T. M., Kjeldsen, H., McCarthy, J. K., Peri, 

M. L., et al., 1993, AJ, 106, 2441. 
Girard, T. M., Wu, H., Lee, J. T., Dyson, S. E., van Altena, W. F., 

et al., 2000, AJ, 119, 2428. 
Guenther, D. B., Kallinger, T., Reegen, P., Weiss, W. W., 

Matthews, J. M., Kuschnig, R., Moffat, A. F. J., Rucinski, 

S. M., Sasselov, D., & Walker, G. A. H., 2007, Commun. 

Asteroseismology, 151, 5. 
Harvey, J., 1985, In Rolfe, E., & Battrick, B., editors, Future 

Missions in Solar, Heliospheric & Space Plasma Physics, ESA 

SP-235, page 199. 
Hatzes, A. P., Guenther, E., Kurster, M., & McArthur, B., 2003, 

In Fridlund, M., Henning, T., & Lacoste, H., editors, Earths: 

DARWIN/TPF and the Search for Extrasolar Terrestrial 

Planets, volume 539, page 441. 
Kambe, E., Ando, H., Sato, B., et al., 2008, PASJ, 60, 45. 
Kervella, P., Thevenin, F., Morel, P., Berthomieu, G., Borde, P., 

& Provost, J., 2004, A&A, 413, 251. 
Kim, K.-M., Han, I., Valyavin, G. G, Plachinda, S., Jang, J. G, 

et al., 2007, PASP, 119, 1052. 
Kjeldsen, H., & Bedding, T. R., 1995, A&A, 293, 87. 
Kjeldse n, H., Bedding, T . R., Arentoft, T., et al., 2008, ApJ. in 

press, arXiv:0804.1182 
Leccia, S., Kjeldsen, H., Bonanno, A., Claudi, R. U., Ventura, R., 

& Paterno, L., 2007, A&A, 464, 1059. 
Marchenko, S. V., 2008, A&A, 479, 845. 

Martic, M., Lebrun, J.-C, Appourchaux, T., & Korzennik, S. G., 

2004, A&A, 418, 295. 
Martic, M., Schmitt, J., Lebrun, J.-C, Barban, C, Connes, P., 

Bouchy, F., Michel, E., Baglin, A., Appourchaux, T., & Bertaux, 

J.-L., 1999, A&A, 351, 993. 
Matthews, J. M., Kuschnig, R., Guenther, D. B., et al., 2004, Nat, 

430, 51. Erratum: 430, 921. 
Mosser, B., Bouchy, F., Catala, C, et al., 2005, A&A, 431, L13. 
Mosser, B., Bouchy, F., Martic, M., Appourchaux, T., Barban, C, 

et al., 2008a, A&A, 478, 197. 
Mosser, B., Deheuvels, S., Michel, E., Thevenin, F., Dupret, M. A., 

Samadi, R., Barban, C, & Goupil, M. J., 2008b, A&A. in press 

(2008arXiv0804.3119). 
Mosser, B., Maillard, J. P., Mekarnia, D., & Gay, J., 1998, A&A, 

340, 457. 

Palle, P. L., Roca Cortes, T., Jimenez, A., Golf, & VIRGO 
Teams, 1999, In Gimenez, A., Guinan, E. F., & Montesinos, B., 
editors, Proc. Workshop on Stellar Structure theory and Tests of 
Convective Energy Transport, volume 173, page 297. ASP Conf. 
Ser. 

Regulo, C, & Roca Cortes, T., 2005, A&A, 444, L5. 

Rupprecht, G., Pepe, F., Mayor, M., et al., 2004, In Moorwood, A. 
F. M., & Masanori, I., editors, Ground-based Instrumentation 
for Astronomy, volume 5492 of Proc. SPIE, page 148. 

Teixeira, T., Kjeldsen, H., Bedding, T. R., et al., 2008, A&A. to 
be submitted. 

Tingley, B., Grundahl, F., & Kjeldsen, K., 2008, MNRAS. to be 
submitted. 

Toutain, T., & Appourchaux, T., 1994, A&A, 289, 649. 



Oscillations in Procyon. I. Observations 



7 



Ulrich, R. K., Garcia, R. A., Robillot, J.-M., Turck-Chieze, S., 
Bertello, L., Charra, J., Dzitko, H., Gabriel, A. H., & Roca 
Cortes, T., 2000, A&A, 364, 799. 



Arentoft et al. 



TABLE 1 
Participating Telescopes 



Identifier 


Telescope /Spectrograph 


Observatory 


Technique 


Ref. 


HARPS 


3.6m/HARPS 


ESO, La Silla, Chile" 


ThAr 


1 


CORALIE 


1.2 m Euler Telescope / CORALIE 


ESO, La Silla, Chile 


ThAr 


2 


McDonald 


2.7m Harlan J. Smith Tel./coude echelle 


McDonald Obs., Texas USA 


iodine 


3 


Lick 


0.6 m CAT/Hamilton echelle 


Lick Obs., California USA 


iodine 


4 


UCLES 


3.9 m AAT/UCLES 


Siding Spring Obs., Australia 


iodine 


4 


Okayama 


1.88m/HIDES 


Okayama Obs., Japan 


iodine 


5 


Tautenburg 


2 m/coude echelle 


Karl Schwarzschild Obs., Germany 


iodine 


6 


SOPHIE 


1.93m/SOPHIE 


Obs. de Haute-Provence, France 


ThAr 


7 


EMILIE 


1.52 m/EMILIE+ AAA 


Obs. de Haute-Provence, France 


white light with iodine 


8 


SARG 


3.58 m TNG/SARG 


ORM, La Palma, Spain 


iodine 


9 


FIES 


2.5 m NOT/FIES 


ORM, La Palma, Spain 


ThAr 


10 



1. Rupprecht et al. (2004); 2. Bouchv & Carrier (20021 ; 3. Endl et al. (20051 : 4. Butler et al. (19961 : 

5. Kambe et al. (2008) ; 6. Hatzes et al. (20031 77. Mosser et al. (200 8a'): 8. B ouchv et al (2002T and paper in prep.; 

9. Claudi et al. (2(555) : 10. FrandseiHT Lindberg (2000) 

"Based on observations collected at the European Southern Observatory, La Silla, Chile (ESO Programme 078.D-0492 (A)). 



TABLE 2 

Summary of observations 



Identifier 


Nights Eff. Obs. 
Allocated 


time 

00 


Spectra 


Median t exp 
(s) 


Deadtime 

00 


/Nyq 

(mHz) 


Noise Level 
(cms" 1 ) 


Noise per Minute 
(ms- 1 ) 


HARPS 


8 


52.0 


5698 


5 


31 


13.8 


2.0 


0.64 


CORALIE 


6 


27.0 


936 


40 


85 


4.0 


9.8 


2.2 


McDonald 


6 


16.2 


1719 


17 


15 


15.6 


14.4 


2.5 


Lick 


14 


95.4 


1900 


180 


54 


2.1 


10.9 


4.7 


UCLES 


12 


41.4 


2451 


16 


45 


8.2 


6.6 


1.9 


Okayama 


20 


83.4 


1997 


54 


110 


3.1 


8.0 


3.2 


Tautenburg 


21 


14.6 


494 


60 


65 


4.0 


22.8 


3.8 


SOPHIE 


9 


35.1 


3924 


23 


28 


9.7 


4.7 


1.2 


EMILIE 


4 


25.7 


1631 


47 


17 


7.8 


10.1 


2.2 


SARG 


4 


15.2 


693 


19 


65 


5.9 


12.6 


2.1 


FIES 


10x1 


12.6 


1087 


18 


45 


7.9 


21.7 


3.4 
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Fig. 1. — Velocity measurements of Procyon, color-coded as follows: HARPS = red; CORALIE = brown; McDonald = gray; Lick = 
cyan; UCLES = blue; Okayama = green; Tautenburg = black; SOPHIE = dark green; SARG = dark blue; FIES = magenta; EMILIE = 
yellow. 

(a) The full time series, before any removal of slow trends (EMILIE and FIES are not shown). 
(6) Close-up of the central ten days (FIES not shown). 

(c) Close-up of a five- hour segment during which three spectrographs observed simultaneously: HARPS (red circles), SOPHIE (dark green 
squares) and SARG (dark blue triangles). All three series have been high-pass filtered to remove slow trends and the SOPHIE and SARG 
data have been smoothed slightly (using a boxcar with a width of three data points). 

(d) The time series of the final noise-optimized uncertainties, showing all 11 telescopes. 
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Fig. 2. — Steps in the adjustments of weights, illustrated using HARPS data from a single night, (a) The velocities, with the slow 
variations removed. (6) The uncertainties, after scaling to satisfy Eq. [T]but before any further adjustments, (c) The uncertainties after 
filtering to remove power on the timescale of the oscillations, (d) The final uncertainties, after adjusting to optimize the noise (see text). 
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Fig. 3. — Same as Fig. [2] but for a single night from EMILIE. 
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Fig. 4. — Time series of velocity measurements of the Sun obtained over 21 days with the GOLF instrument on the SOHO spacecraft. 
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Fig. 5. — Power spectra for all 11 telescopes, together with that of the combined series. Note that the vertical scale is not the same for 
all panels. 
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Frequency (mHz) 



Fig. 6. — Final power spectrum based on the noise-optimized weights (lower panel), and also without applying the weights (upper panel). 
The inset shows the spectral window. 




10 100 1000 10000 

Frequency (/CiHz) 

Fig. 7. — Power density spectrum of Procyon from the HARPS data, and the same after smoothing. The lower two dashed lines show 
the solar activity level at minimum and maximum, and the upper line is the solar maximum activity shifted upwards by a factor of 40. 
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Frequency (/iHz) 

Fig. 8. — Smoothed power density spectra (see FigO for the nine telescopes shown in Fig.[l]i, showing a similar background from stellar 
noise at low frequencies and different levels of white noise at high frequencies. The color coding is the same as in Fig.[T] 
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Fig. 10. — Smoothed amplitude curves for Procyon from ten 2-day segments of the combined time series (thin lines), together with their 
mean (thick line). 




